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Abstract — In the landmark paper 1 1 1 by Hassibi and Hochwald, 
it is claimed without proof that the upper triangular matrix R 
encountered during the sphere decoding of any linear dispersion 
code is full-ranked whenever the rate of the code is less than the 
minimum of the number of transmit and receive antennas. In this 
paper, we show that this claim is true only when the number of 
receive antennas is at least as much as the number of transmit 
antennas. We also show that all known families of high rate 
(rate greater than 1 complex symbol per channel use) multigroup 
ML decodable codes have rank-deficient R matrix even when 
the criterion on rate is satisfied, and that this rank-deficiency 
problem arises only in asymmetric MIMO with number of 
receive antennas less than the number of transmit antennas. 
Unlike the codes with full-rank R matrix, the average sphere 
decoding complexity of the STBCs whose R matrix is rank- 
deficient is polynomial in the constellation size. We derive the 
sphere decoding complexity of most of the known high rate 
multigroup ML decodable codes and show that for each code, 
the complexity is a decreasing function of the number of receive 
antennas. 

I. Introduction 

A. System Model and Definitions 

We consider Space-Time Block Codes (STBCs) for an 
N transmit antenna, M receive antenna, quasi-static MIMO 
channel (N x M MIMO system) with Rayleigh flat fading. 
The system can be modeled as 



Y = XH + W, 



(1) 



where X is the T x N codeword matrix transmitted over T 
channel uses, Y is the T x M received matrix, H is the 
N x M channel matrix and the T x M matrix W is the 
additive noise at the receiver. The entries of H and W are 
i.i.d. zero mean, unit variance, circularly symmetric complex 
Gaussian random variables. 

Definition 1: (STBC) An STBC C encoding K real inde- 
pendent information symbols, denoted by xi, i = 1, • • ■ , K, 
is a set of complex matrices given by 



C 




(2) 



where the T x N complex matrices Ai, . . . ,Ak, which 
are called linear dispersion or weight matrices, are linearly 
independent over the real field R ffl. Ell, and the finite set 
A C M. K is called the signal set. 

Definition 2: (Code Rate) The rate of an STBC is the 
average number of information symbols transmitted in each 



channel use. For the STBC given by (0, the code rate is 
R = real symbols per channel use, or R = complex 
symbols per channel use (cspcu). 

The linear independence of the weight matrices in the 
definition of an STBC implies that R < N. Throughout this 
paper, unless otherwise specified, the code rate is taken to be 
in terms of complex symbols per channel use. 

Generally, the signal set A is chosen in such a way that the 
STBC C has full-diversity and large coding gain. In most cases 
A is chosen to be a subset of 0Z , where 6 R KxK is 
a full-ranked matrix. One such instance is when the symbols 
are partitioned into multiple encoding groups, and each group 
of symbols is encoded independently of other groups using 
a lattice constellation, such as in Clifford Unitary Weight 
Designs Q (in which case is an orthogonal matrix), Quasi- 
Orthogonal STBCs [3 1 and Coordinate Interleaved Orthogonal 
Designs 0. There are also instances where the real symbols 
are encoded independently using regular PAM constellations 
of possibly different minimum distances 0, in which case 
is a diagonal matrix with positive entries. 

For a complex matrix A, let its real and imaginary com- 
ponents be denoted by A] and Aq, respectively. Let vec(A) 
denote the complex vector obtained by stacking the columns 
of A one below the other and 

vec(A) = [uec(A/) T vec(AQ) T ] T . 

Now, the system model given by (Q~|) can be expressed as 

y = vec(Y) — Gx + w, 

where x = [xi, . . . , xk] t , W = vec(W) and the equivalent 
channel matrix G € M. 2MTxK is given by 



[uec(A]H) «ec(A 2 H) 



"(AkH)]. 



Consider the vector of transformed information symbols 
s = _1 x which takes values from A' = ~ 1 A C Z , where 
Z denotes the ring of integers. The components of s take finite 
integer values, i.e., Sj <E A q C Z, with \ A q \ = q for some finite 
q. Hence one can use a sphere decoder ||6] to decode s and 
then obtain the ML estimate of the information vector x. The 
ML decoder output is given by 



where 



s = arg min lly — G0s|| w , 
seA' 

\f denotes the Frobenius norm of a matrix. 



(3) 



B. Motivation for our results 



, Ai 7 of the STBC in fljj) are as follows: 



It is claimed in [1| without proof that R < min{M, N} is 
a sufficient condition for the system of equations defined by 
d3) to be not underdetermined, i.e., for 

rank{G) = rank{G@) = K, with K < 2MT. 

In Section [ill we show that the claim made in [T) is true only 
for M > N. This observation is the gateway to the new results 
presented from Section HT1 onwards. 

For a system where rank(G) — K, the sphere decoder 
complexity, averaged over noise and channel realizations, 
is independent of the constellation size q and is roughly 
polynomial in the dimension of the sphere decoding search 
GD, 0, 0. However, if the rank of G is less than K, the 
average sphere decoding complexity is no more independent 
of the constellation size. When rank(G) = K' < K, the con- 
ventional sphere decoder needs to be modified as follows ]9] : 
The R matrix resulting from the QR-decomposition of G© 
has the form R = [R a R 6 ] e R K ' xK , where R a is a K' x K' 
upper triangular full-rank matrix. There is a corresponding 
partition of s as [s„ s^] T . If q is the size of the regular 
PAM constellation used, then for each of the q K ~ K values 
of Sfc, the conditionally optimal estimate of s a can be found 
by first removing the interference from %, and then using a 
sphere decoder with the upper triangular matrix R a to obtain 
an estimate of s . Then, from the resulting q K ~ K estimates 
of s, the optimal vector is chosen. This observation leads to 
the following lemma. 

Lemma 1: For an STBC whose equivalent channel matrix 
G is such that rank(G) = K' < K, the sphere decoding 
complexity as a function of q is of the order of q K ~ ' '• ■• 
O 



In the rest of the paper, by sphere decoding complexity we 
mean the average sphere decoding complexity and the focus is 
on the dependence of the complexity on the constellation size 
q and not on the dimension of the sphere decoding search. We 
now introduce the notion of singularity of an STBC which is 
a direct indicator of its sphere decoding complexity. 

Definition 3: Let C be an STBC with rate R = £p and let 
M > R. We say that C is singular for M receive antennas if 
rank(G) < K with probability 1. Otherwise, C is said to be 
non-singular for M receive antennas. 

To illustrate the effect of the rank of G and singularity of 
an STBC on the sphere decoding complexity as a function of 
the constellation size, consider the following example. 

Example 1: In 1 10], a rate R = code for N = 4 antennas 
with K = 17 real symbols was presented. Define X 



and Z 



1 
[l 

Kronecker product by A 



1 

, and for any two matrices A, B, denote their 
B. Then, the 17 weight matrices 



Ai = 


I2OI2, 


A 2 = 


= iL(g 


I 2 , 


A 3 = 


ZX®I 2 , 


A 4 = 


tX®Z, 


A 5 = 


= X(g> 


ZX, 


A 6 - 


iX®X, 


A 7 = 


I 2 ® zx, 


A 8 = 


= iZ(g 


x, 


A 9 = 


ZX®X, 


Aio = 


= iZX ® ZX, 


An 


= zx 


®x, 


A12 = 


= iX <g) I 2 


A13 = 


= Z <g> ZX, 


A14 


= ih 


ox, 


A15 = 


= iZ <S) Z, 


Ale = 


= H2 ® I2 and 


A17 


= 0.2 


<g>Z. 







Let the number of receive antennas be M = 3. Then, M > R 
and the equivalent channel matrix G is of size 24 x 17. Now 
consider the following randomly generated channel matrix 



0.3457 + 0.2299i 
0.7316 - 0.5338i 
0.5140 + 0.9689i 
-0.2146 - 1.2102i 



0.2078 
-0.5567- 

0.6282- 
-0.8111 - 



0.0723i 
0.1707i 
0.2257i 
0.2212i 



-0.7558 - 0.6116i 
-0.5724 -0.0212i 
-2.0819 -0.1166i 
1.0171 + 0.4439i 



The resulting G matrix has rank only 16 and the structure 
of the upper triangular matrix R of size 24 x 17 obtained 
upon its QR-decomposition is shown in (O at the top of 
next page. The non-zero entries of R are denoted by 'a'. It is 
clear that the first 16 columns of R are linearly independent 
and the last column lies in the span of the first 16 columns. 
Removing the last 8 rows of R, which are all zeros, we 
get the 16 x 17 real matrix R' which is used by the sphere 
decoder to find the ML estimate of the information vector. In 
this case, from Lemma Q] the complexity of sphere decoding 
this STBC for this particular channel realization H is of the 
order of q 17 ~ le = q. In Section IIII-AI by establishing that 
rank(G) = 16 with probability 1, we show that this STBC 
is singular for M = 3. 

C. Contributions of the Paper 

The contributions and organization of this paper are as 
follows. 

• We introduce the notion of singularity of an STBC which 
is a direct indicator of its sphere decoding complexity. 

• We show that contrary to the claim made in (T), 
R < min{M, N} is not a sufficient condition for an 
STBC to be non-singular for M receive antennas (The- 
orem Q] Section [El). We show that the case of singular 
STBCs arises only for asymmetric MIMO systems with 
M < N (Proposition Q] Section [TTJ. 

• We show that all known families of high rate {R > 1) 
multigroup ML decodablef] codes HO), fl~Q, Q3, C21 
and [ 14] are non-singular for certain values of M {Sec- 
tion|lIIl see Table U for a summary of results). 

« We derive the sphere decoding complexity of almost all 
known high-rate multigroup ML decodable codes and 
show that in each case the sphere decoding complexity is 
a decreasing function of the number of receive antennas 
M (Section [nil see Table [J). We show that even when 
an STBC is singular, multigroup ML decodability helps 
reduce the sphere decoding complexity. The reduction in 

'An STBC is g-group or multigroup ML decodable if its symbols can be 
partitioned into g groups and each group of symbols can be ML decoded 
independent of others. 
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(4) 



complexity is from O 

\Q K ~ K ) to O (q^—y where g 
is the number of ML decoding groups (Section HlB . 
Some related open problems are discussed in Section [IV] 
Notations: Throughout the paper, matrices (vectors) are 
denoted in bold, uppercase (lowercase) letters. For a complex 
matrix A, its transpose, conjugate and conjugate-transpose 
are denoted by A T , A and A H , respectively. For a square 
matrix A, det(A) denotes its determinant. The n x n identity 
matrix is denoted by I„ and is the null matrix of appro- 
priate dimension. Unless used as a subscript or to denote 
indices, i = y/—X. For square matrices Aj, j = l,...,d, 
diag (Ai, Ad) denotes the square, block-diagonal matrix 
with Ai, . . . , Ad on the diagonal, in that order. The field of 
complex numbers is denoted by C. 

II. Basic results on the rank of the equivalent 

CHANNEL MATRIX 

We present a few results which we will use in the following 
sections to derive the rank of G for specific STBCs. The 
following result shows that if any STBC is singular for M 
receive antennas, then M < N. Thus, the rank-deficiency 
problem arises only in asymmetric MIMO with M < N. 

Proposition 1: If M > N, eveiy T x N STBC is non- 
singular for AI receive antennas. 

Proof: Since vec(-) is an isomorphism from the IE- 
vector space C TxM to R 2MT , it is enough to show that 
AiH, . . . ,AifH are linearly independent with probability 1. 
Suppose H is full-ranked, i.e., rank(H) — N, then there 
exists a matrix H f G C MxN such that HH f = l N . If V = 



2fc=i a «AiH = 0, it would mean that VH^ = J2f=i a i^-i = 0. 
Since Aj are the weight matrices of an STBC, they are 
linearly independent and hence a, = 0, i = 1, . . . , K. Thus 
ranfc(G) = K if H is full-ranked. Since H is full-ranked with 
probability 1 [15], we have shown that any STBC is non- 
singular for M receive antennas if M > N. ■ 

Let (Ai, A2, . . . , Ak) denote the R-linear subspace of 
C TxN spanned by the matrices Ai, . . . ,Ak- 

Proposition 2: Let Bj., . . . , Br be T x N complex matrices 
such that (Ai, . . . ,A K ) = (Bi, . . . ,B X > and let H € C NxM 
be any matrix. Then the column spaces of the matrices 

G A (H) = [«ec(A]H) uec(A 2 H) • • • vec(A K B)} and (5) 
G B (H) = [uec(B x H) wc(B 2 H) ■ • • vec(B K B)] (6) 

are identical. In particular, ranfc(GA(H)) = ranfc(GB(H)). 

Proof: Let v be any vector in the column space of Ga(H). 
Then, v = 2~2^=i a-iwec(AiH), for some choice of real numbers 
ai, i = 1,...,K. Since each of the Aj € (Bi, . . . ,B#), 
every Aj can be written as some real linear combination of 
Bi,...,Bjf. It follows that every uec(AjH) can be written 
as some real linear combination of uec(BiH), . . . , uec(B^-H). 
Hence v belongs to the column space of Gb(H). Similarly we 
can show that every vector in the column space of Gb(H) 
belongs to the column space of Ga(H) also. This completes 
the proof. ■ 
The following result shows that if every weight matrix of a 
given STBC is multiplied on the left by a constant invertible 
matrix, then the rank of the equivalent channel matrix is 
unchanged. 



Proposition 3: Let C € C TxT be any full-rank matrix, 
Ri = CAj, i = 1, . . . , K and H be any NxM complex matrix. 
Then we have rank (Ga(H)) = rank (Gb(H)), where Ga(H) 
and Gb(H) are as defined in (0 and ©. 

Proof: It suffices to show that the subspaces 
(A]H, . . . , A^-H) and (CAiH, . . . , CA^-H) have the same 
dimension. The proof is complete if we show that the vector 
space homomorphism (p from the former subspace to the latter, 
that sends V = 2~2i=i a iAiH to Yli=i QtCAjH is a one to one 
map. If V e ker((p), then ip(V) = C J2?=i a*A;H = CV = 0. 
Since C is invertible, this means that V = 0. This completes 
the proof. ■ 

It is shown in [14] that for any N > 1, there exists 
an explicitly constructable set of N 2 matrices belonging to 
C NxN that are unitary, Hermitian and linearly independent 
over R. This set of matrices forms a basis for the space of 
N x N Hermitian matrices. Denote by C§ erm any STBC 
obtained by using these N 2 matrices as weight matrices. For 
positive integers n and m, define the function 

f(n, m) = n 2 — ((n — m) + ) 2 , 

where (a) + = max{a,0}. We now state the main result of 
this paper. 

Theorem 1: The equivalent channel matrix of the STBC 
QHerm f or j^j rece j ve antennas has rank f(N, M) with 
probability 1. 

Proof: See Appendix A. ■ 
The rate of C^ erm is y and the rank of the equivalent 
channel is less than K = N 2 whenever M < N. Thus, this 
STBC is singular for all f < M < N. 

Example 2: Consider the STBC C^ erm used in an asym- 
metric MIMO system with AI = 2 receive antennas. In this 
case, i?,= |< min{M, N} and the equivalent channel matrix 
G is of size 12 x 9. From Theorem QJ we know that the rank 
of G is equal to /(3, 2) = 8 with probability 1. Hence, this 
STBC is singular for 2 receive antennas. The 9 weight matrices 
of the STBC Cf erm are as follows 
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1 








, A 5 = 
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The structure of the 12 x 9 upper triangular matrix R' 
obtained from the QR decomposition of G when H equals 
the following randomly generated matrix 

"-0.5688 - 0.8117* -0.1723 + 1.8282i" 
0.4926 + 0.0742i 0.1525 - 0.4716i , 
0.5905 + 0.5107i -0.8244 + 0.1325i 



is given by 
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The sphere decoder uses the 8x9 matrix R obtained from 
R 1 by deleting its last 4 rows which are all zero. Hence, 
for this particular channel realization H, the sphere decoding 
complexity is of the order of q 9 ~ 8 = q. 

The remaining part of this section is concerned with multi- 
group ML decodable codes. Suppose the code obtained from 
an STBC C with a signal set A is g-group ML decodable for 
some g > 1. The information symbols {x%, . . . , xk} can be 
partitioned into g vectors xj l; . . . ,xi of length Ai, . . . , X g 
respectively such that each symbol vector can be ML decoded 
independently of other symbol vectors. There is a corre- 
sponding partition of the channel matrix into g submatrices 
Gi, . . . , G ff , such that 

9 

Gx = ^G fc x Ifc , G fc eR 2MTxAfc , for k = l,...,g. (7) 

k=l 

In Theorem 2 of lfl6l it is shown that for any k ^ k 1 and 
any channel realization H, every column of Gfc is orthogonal 
to every column of Gfc<. As a direct consequence of this, we 
have the following proposition. 

Proposition 4: For any g-group ML decodable STBC and 
any channel realization H, rank(G) — 2~2t=i rank{Gk)- 

Proof: Since the column spaces of G^, k = l,...,g are 
orthogonal to each other, the dimension of the column space 
of G is equal to the sum of the dimensions of the column 
spaces of G&, k = 1, . . . , g. ■ 

III. Sphere decoding complexity of some known 

FAMILIES OF CODES 

In this section, we show that all known families of high- 
rate (R > 1) multigroup ML decodable codes are singular for 
certain number of receive antennas. Using the properties of the 
rank of the equivalent channel matrix derived in the previous 
section, we now derive the sphere decoding complexities 
of these known multigroup ML decodable STBCs. Table U 
summarizes the results of this section. The table lists the 
sphere decoding complexity and the minimum number of 
receive antennas for non-singularity of C^ erm , the codes 
in iflOl . Ifl2l . fl3l . and the codes in lfl4l corresponding to 
even number of ML decoding groups. 



TABLE I 

Comparison of Sphere Decoding Complexities of known singular STBCs 



Code 


Transmit 
Antennas 

AT 
JV 


Delay 
l 


Groups 
9 


Rate 
(cspcu) 

H 


Minimum M 
for 

non-singularity 


Order of Sphere 
Decoding Complexity 
Exponent of q''* 


C He.rm ( Theorem [JJ 


> 1 


N 


1 


V 
2 


JV 


((N-M)+) 2 


Ren et. al. QO) 


4 


4 


2 


17 

T 


4 






((4-M)+) 


2 






Ren et al. CD 


> 1 


> 2jV 


2 


JV ^T 1 


jV 


(N - M)+ ■ (T — N — M) 


Srinath et al. 1 13 ] 


2 m , m > 2 


N 


2 


V , 1 

4 JV 


\' 

2 






2 


Natarajan et. al. |14|t 


1 9 — 1 1 

ns2 L ^~ J , n > 1 


N 


2£, e > 1 


JV . B 2 -S 


JV 


( 


I » 








2 

) 


g 2 3-l ^ 2JV 






2 - 


n2 L - £ 2 J "J, n > 1 


gN 


2£, £ > 1 


JV i 9-1 


JV 


( 




t- 2^J 




+ ) 


2 


2 9-l 2JV 


29-^ 




-J 



t The size of the real constellation used is denoted by q. 

* For any real number a, a + is defined as max{a, 0}. 

t 1 1 4 1 contains codes for all g > 1 and not just even values of g. 



A. Fast-group-decodable STB C from Ren et. al. MOV 

In lITOl . a 2-group decodable STBC for 4 transmit antennas 
with R = 17/8 was constructed. For this code, with the nota- 
tions as used in (01, Xx 1 = Xi, Ai = 1, xj 2 = [x2, . . . , xiq] t , 
A2 = 16 and G = [Gi G2]. The weight matrix corre- 
sponding to x\ is I4. Since any two weight matrices from 
different groups are Hurwitz-Radon orthogonal, i.e., satisfy 
hf Aj + Ai = 0, all the matrices in the second group must 
be skew-Hermitian, i.e., A^ = —A4, i = 2, ■ ■ • , 17. As a 
result, we can express these weight matrices as Aj = iLiA^, 
i = 2, • • • , 17, with A^ being Hermitian matrices. It is clear 
that for any arbitrary channel realization H, Gi = [wec(H)] 
which is non-zero with probability 1. Hence, rank(G\) = 1 
with probability 1. From Proposition [3] and Theorem [TJ 
rank(G 2 ) = /(4,M) with probability 1. So, 

1) for AI — 3, rankiGz) = 15. Since xx 1 and xj 2 can be 
decoded independently of each other, from Proposition 
[5] the sphere decoding complexity of the first group 
is independent of q while that of the second group is 
0(g( 16 ~ 15 )) = 0(q). Consequently, the sphere decoding 
complexity of the code in ifTOl is 0(q) for 3 receive 
antennas. 

2) For M > 4, /(4, M) = 16 and the STBC in d is 
non-singular for 4 or more receive antennas. Hence, its 
sphere decoding complexity is independent of q. 

B. High-rate 2-group ML decodable codes from Srinath et. 

ai. an} 

A family of 2-group ML decodable STBCs was constructed 
in US for N = 2 m , m > 1 antennas with rate R = f + i 
cspcu. This family includes the rate | code of ifTTI for N = 4 
as a special case. The number of symbols in the STBC is 
K = ^- + 2. In the rest of this subsection we show that the 

'M-M) + y 



sphere decoding complexity is O [ q 



which is 



\a\ is the smallest integer greater than or equal to a. Note 
that the sphere decoding complexity is a decreasing function 
of the number of receive antennas M. 

Derivation of sphere decoding complexity: The STBCs 
constructed in 1131 have a block diagonal structure. The weight 
matrices for N = 2 m antennas are of the form diag(Vi, V2), 
where Vi, V2 € C" x ™ are unitary and n = Noting that the 
STBC is 2-group decodable, denote the set of weight matrices 
belonging to the first and the second group by Q\ and Q2, 
respectively. For all the matrices in Q\, Vi is constant, say 
Fi, and for all the matrices in Q2, V2 is constant, say F2. 
Each group contains n 2 + 1 real symbols. We will now derive 
the rank of the submatrix Gi of G that corresponds to Q\. 

Let Q\ = {Ai, . . . , A„2 +1 }. Consider the set 

G[ ={A' l =CA 4 ,z = l,...,n 2 + l}, 

where C = diag(l n , $2 )• Since C is unitary, from Proposi- 
tion [3] the rank of Gi equals the rank of G' l5 the equivalent 
channel matrix corresponding to Q[. Since the multiplication 
of all the weight matrices of an STBC by a unitary matrix 
does not affect its multigroup ML decodability, for any matrix 

+ 1, we have 

0, (8) 



B belonging to Q2 and any j = 1, . . . , r 
Af (CB) + (CB) ff A; 
where A' and CB are block diagonal and are of the form 



A' 



Fi 





Do 



and CB 







for some unitary matrices D 3 and V. Since A'j and CB 
satisfy ([H} and are block diagonal, from Lemma 1 of [13|, 



Dfa„ 



(il n ) ~Dj = 0. Thus, for j = 1,. 



1, 



Dj is a 



polynomial in q, and so is large for all \R] < M < f , where 



n x n unitary, Hermitian matrix. Next we use Proposition |2] 
to find the rank of G^. 

Note that (A' 1; . . . , A^ 2+1 ) is same as the span of 



(9) 
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D„2+l_ 



Since Ai, . . . , A n 2 +1 are linearly independent, A' l5 . . . ,A' n2+1 
are also linearly independent. Further, the first matrix in (0 
is linearly independent of the remaining matrices and hence 
the dimension of the span of the remaining matrices in (0 is 
n 2 . Without loss of generality, let us assume that Di, . . . , D n 2 



are linearly independent, thus (Di, . . . , D„2 
all n x n Hermitian matrices. Then, (A^, . . 
the space spanned by 



"Fi 










" 
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D„2_ 



is the space of 
, A^ 2+1 ) equals 



(10) 



From Proposition [2] it is enough if we concentrate on 
the STBC whose weight matrices are given by ([ToV Let 

"Hi 



the channel matrix be partitioned as H = 

ixM 



H 2 



where 



Hi,H2 G C™ . We need to compute the dimension of the 
space spanned by the weight matrices multiplied on the right 
by H which is 



FiHi" 

















DH 2. 




D„2H 2 



With probability 1, Hi is non-zero and hence the first ma- 
trix is linearly independent of the remaining matrices. From 
Theorem Q] the dimension of the span of the remaining n 2 
matrices is f(n, M) = n 2 — {{n — M ) + ) 2 with probability 
1. Thus ranfc(Gi) equals f(n,M) + 1 with probability 
1. A similar result can also proved for the second ML 
decoding group, i.e., for ranfc(G 2 ). From Proposition |4] 
rank(G) = rank (Gi) + ranfc(G 2 ) which equals 



K' 



Comparing this with K = 2 L^- + lV we see that the STBC 
given in |fl3ll is non-singular only if M > % ps 2R. Hence the 




code is singular for all \R\ < M < y. 

Now, the two groups of symbols can be ML decoded 
independently of each other, and the number of symbols in 
each group is -y, with rank{G\) — rank{G\) = Thus, 

the sphere decoding complexity of the STBC is O (q ~ 

instead of O (q K ~ K ^. Hence, multigroup ML decodability 
reduces the sphere decoding complexity even if the STBC is 
singular. 

C. Two group ML decodable codes from Ren et. al. [12 / 

In fl2l . 2-group ML decodable codes for all N > 1 and all 
even T > 2N were constructed with rate R~N — N ~ x . 
The number of symbols per group is = TN — N 2 + 1. 
In this subsection, we show that the codes of this family 
are singular for all \R~\ < M < N receive antennas and that 
their sphere decoding complexity is O ^g( Ar - M ) + -( T - A '- M )^ , 
Using Proposition Q] it is clear that these codes are non- 
singular if and only if M > N. 



Derivation of sphere decoding complexity: The structure 
and derivation of the sphere decoding complexity of these 
codes is similar to that of the codes from iFPTl . which was 
discussed in Section ITlI-BI The weight matrices of the STBCs 



in |[T2| are of the form 



Vi 

v 2 



-xN 



For 



where Vi,V2 € 

all the matrices in the firsf group, Vi = ±Fj for some 
constant matrix Fj, and for all the matrices in the second 
group, V2 = ±F 2 for some constant matrix F 2 . The STBCs 
constructed in lfl2l are such that for each k — 1,2, the Vfc 
submatrices of any two weight matrices belonging to different 
groups are Hurwitz-Radon orthogonal. We consider the case 
where F^, k = 1,2 are semi-unitary i.e., F^F^ = Ijy. We 
derive the sphere decoding complexity only for the first group. 
Using a similar argument the complexity for the second group 
can be derived, and it is same as that of the first group. 
Since F 2 is semi-unitary, there exists a unitary T x T matrix 



F~ 2 such that F 2 F 2 







. Consider the new STBC C 



obtained by multiplying all the weight matrices of the original 





STBC on the left by C = 



It 




Then, the lower 



submatrix of all the weight matrices of the second group are 
~* A Since the lower submatrix of every matrix 



of the form 







in the first group is Hurwitz-Radon orthogonal to 



±Ijv 




, the 



weight matrices in the first group of C have the following 
fVil 

, where Vi = ±Fi, B is an N x N Hermitian 



structure: 



matrix and E 6 ~ NxN . Let Bi, . . . , B^ 2 be an y basis for 
the space of N x N Hermitian matrices over R, L = TN—2N 2 
and Ei, . . . ,Ei be 



1 







i 







1 















i 































which is a standard basis for the space of (f - N) x N 
complex matrices over R. Clearly the space spanned by the 
weight matrices of the first group of C is a subspace of 
the space spanned by the following linearly independent 



matrices: 






"Fi" 




" 









B„ 













forn = l,...,7V 2 



for I = 1, . . . , L. 

(11) 



From dimension count, it is clear that the matri- 
ces in ( fTTT i form a basis for the space spanned by 
the weight matrices of the first group of C, For 



any non-zero channel realization H G 



spaces Vi 



and V3 




V 2 = 




B„H 




iNxM 



n = 1, 



the sub- R = 



,N 2 



) are such that their pair- 



wise intersections contain only the all zero matrix. Thus, the 
rank of the equivalent channel matrix of the first group of C' 
is 

rank(G 1 ) = dim(Vi) + dim(y2) + dim(Vs). 

With probability 1, dim(Vi) = 1 and from Theorem Q] 
dim(V 2 ) = /(iV, M). It is straightforward to show that 
dim{Vz) = (T - 2N) ■ min{N, M} with probability 1. From 
Proposition Q] we know that every STBC is non-singular for 
M > N. We thus consider only the case M < N. Then from 
Proposition [3] the rank of the equivalent channel matrix of the 
first ML decoding group of the STBC given in lfl2l is 



ran 



fc(Gi) = ranfc(Gi) = (T — 2N) ■ M + f(N, M) + 1. 



Compare this 

K 
2 



of symbols in the first 
1. Thus, for any M < N, 



with the number 
group f = (T - 2N) ■ N + N 2 + 
rank(Gi) < -j with probability 1. 

Using a similar argument it can be shown that the rank of 
the equivalent channel matrix of the second group also equals 
(T - 2N) ■ M + f(N, M) + 1 with probability 1. Since, the 
two groups can be ML decoded independently of each other, 
the complexity of sphere decoding the STBC proposed in ifTZl 

is 0( g (^-M>(T-JV-Af)) for any M < N and the STBC is 

singular for all \R\ < M < N. 

Example 3: Consider the STBC from ifTZl for N = 4 and 
T = 2N = 8. The rate of this code is R = ^ and the number 
of symbols per decoding group is 4- = 17. From the above 
discussion it is clear that this STBC is singular for M = 3. The 
rank of the equivalent channel matrix of each ML decoding 
group equals 16 with probability 1 and hence the sphere 
decoding complexity is 0(q). 

It is interesting to compare this code with the code from iTTOI 
which we have discussed in Section IIII-AI and Example Q] 
Both codes have the same parameters N, R and both have a 
sphere decoding complexity that is linear in the constellation 
size q. However, the code in [ 1 1 is fast-group-decodable and 
5 levels can be removed from the sphere decoding search tree 
of the second decoding group. Hence, after conditioning on 
the value taken by one of the real symbols (to account for 
the reduced rank of the equivalent channel matrix), the code 
in ifTUll uses a 10-dimensional sphere decoder to decode the 
second group of 16 symbols. For decoding each ML decoding 
group, the code from lfT2l . uses a 16 dimensional search tree 
after conditioning on the value of one of the real symbols. 
Thus, the sphere decoding complexity of the code from [ 1 1 
is less than that of [12|. 

D. Multigroup ML decodable codes from Natarajan et. al i!4\I 
In fl4l . delay optimal g-group ML decodable codes were 



g2 ^-i + § 2N 9 ■ ^ n tn ' s subsection we show that the sphere 
decoding complexity of the codes with even g is of the order of 

q V V / / and that the codes are non-singular only 

for M > Also in El, non-delay optimal codes with 

T = gN, N = n2LS-J, n > 1 were constructed. We show 
that the sphere decoding complexity of these codes for even 

((n-2^1 M Y) 2 

values of g is of the order of q V \ / / . Simulation 

results show that the STBCs in [14] for odd values of g are 
also singular for certain number of receive antennas. 

Delay-optimal codes: The number of symbols per group 



is — = ir + g — 1 and let m ~ 2L £ 2 _ J. We will now de- 



K_ 

<j 



rive the rank of the equivalent channel matrix Gi of the 
STBC corresponding to the first group. The weight ma- 
trices of the first group have a block diagonal structure 
diag (Di,D 2 , . . . ,D 9 ), where each e c nrnxnm . The first 
block Di is one of the n 2 matrices of the form V®Ui, where 
V € C" x " i s Hermitian, and the remaining g — 1 matrices 
Dj, j — 2, ■ ■ ■ ,g are of the form ±I„ ® Uj for some set 
of g unitary mxm matrices Ui,...,U 9 . Let us multiply 
all the weight matrices of the first group on the right by 
C = diag (l„ ® Uf , ■ ■ • ,I n ® U^). Clearly, the new set of 
weight matrices A' 1 ,.,.,Ak> also have a block diagonal 

g 

structure diag (Di, D 2 , . . . , D g ) where Di is one of the n 2 
matrices of the form V €5 1„, where V is n x n Hermitian and 
the remaining g — 1 blocks are of the form ±I„ <g> I m . It is 
straightforward to show that there exists a permutation matrix 
P such that P • (V ® I m ) • P T = I m <g> V for any V € C nxn . 

11 11 iTP T*C 

Consider the matrices A, = C AX ,7 = 1,...,—, where 
C' = diag(P,I nm ,...,I nm ). Let Bi,...,B n a be any basis 
for the space of n x n Hermitian matrices. From dimension 
count and the structure of the A" matrices, (Aj,...,Ak) 

equals the span of the matrices Ai, . . . , Ak which given by 

g 

diag(l m <S)Bi,0, . . . , 0), I = 1, . . . , n 2 , diag(0, I nm , 0, . . . , 0), 
. . . , diag(0, . . . , 0, I nm ). Since C /T is unitary, the statistics of 
H and C /T H are same. Along with Propositions [2] and |3] it is 
thus clear that rank (Gi) has the same statistics as the rank of 
the equivalent channel matrix Gi of the STBC whose weight 
matrices are Ai, . . . ,Ak, 

i T 



Let the channel realization be H = [frf 



H 



where Hfc € C 



nrn X M 



T 3 \ 



k = 1, . . . ,g. We are interested in the 



dimension of (A]H, . . . , AkH) which is the span of 







for I = 1, 



" " 




" " 


H 2 












H 9 . 



constructed for all g > 1, N — ng2^ 



n > 1, with rate 



(12) 



Thus, with probability 1, ranfc(Gi) equals the sum of g—1 and 
the dimension of the span of the first n 2 matrices in (fT2l . Let 
us rewrite Hi as [H^ . . . H^ m ] T , where Hij G C nxM 
for j — 1, . . . , m. Thus, the dimension of the span of first n 2 



matrices in (TTZt is same as that of 



given in dl~3l > at the top of the next page. The matrix R' is of the 



B1H1 i 




B2H11 




B n 2Hi,i 


BiHi^ 




B2Hi ; 2 


j ■ ■ • ? 


B„2Hi : 2 


_BiHi jm _ 




_B2Hi im _ 




_B n 2Hi^ m ^ 



This in turn, is equal to the dimension of the span of the 
following matrices 



[BiHi,i BiHi,2 
[B2Hi i i B2Hi i 2 



BiHi^ m ] 
B2Hi /m ] 



[B„2Hi ! i B„2Hi^ 2 ■ • • B„2Hi, m ] . 

From Theorem Q] this dimension equals f(n, mM) with 
probability 1. Hence, 

ranfc(Gi) = rank{G\) = f(n, mM) + g — 1 



with probability 1. Compare this with the number of sym- 
bols per group y = n 2 + g — 1. Thus, the delay optimal 
STBCs of 1 14 1 for even values of g are singular whenever 
\R\ < M < — j^t— ■ Similar results on the rank of the equiv- 

2 L_ 2 _ J 

alent channel matrix can be proved for other ML decoding 
groups also. Thus, the sphere decoding complexity of these 

codes is of the order of J J . This STBC 

is singular only for M > g2 ^_ 2 « 2R. For g = 2, and equal 
values of N, the delay-optimal codes of [14] and the codes 
of lfl3l have equal rate and the same order of sphere decoding 
complexity. 

Example 4: Consider the g = 2, N = 6 delay-optimal code 
of lfl4l . The underlying STBC has rate R = | and there are 
Y = 10 symbols per ML decoding group. From the ongoing 
discussion, for M = 2 receive antennas rank(Gt) = 9 with 
probability 1 for k = 1,2. Hence, this STBC is singular for 
M = 2 and the sphere decoding complexity is 0(q). Let us 
denote the weight matrices of C^ erm given in Example |2] by 
A l5 . . . ,A 9 . Then, the 20 weight matrices of the N = 6, g = 2 
code of [14 1 are given by 

I, ^- 1 ---» 9 » A io- 



A, = 



M+10 



I 3 






iA'„ 



= 1,...,9, and A 20 = 





Is 

-Is 








Now consider the 24 x 20 G matrix corresponding to the 
randomly generated channel realization 



H 



-0.0583+ 1.2105z 
-1.3669- 0.1373i 
-0.3104- 1.5120i 
-1.2690 - 0.5937i 
0.5942 + 0.9578* 
-0.6279 - 0.758H 



0.0708 + 0.6795i 
-0.3850 + 0.0877i 

0.2146+ 1.0159i 
-0.4245- 1.3866i 

0.3465 - 0.1398* 

0.5228 -0.8541* 



The rank of G is only 18 and the structure of the 24 x 20 upper 
triangular matrix R' obtained from the QR decomposition is 



form 



Ri 




o 9 > 


;10 


0i> 


:10 


0i> 


;10 


o 9 > 


;10 


R 2 




o 5 > 


;10 


o 5 > 


;10 



where Ri,R2 G 



The 10 x 10 



all zero submatrix at the upper right corner of R' is due to the 
2-group ML decodability property of the code. The two sphere 
decoders corresponding to the two ML decoding groups use 
the matrices Ri and R2 respectively. Clearly the rank of Ri 
and R 2 is 9, and hence the sphere decoding complexity is 
0(q). 

Non delay-optimal codes: The non delay-optimal code 
for N = 7x2 L^ - J has rate R — + and number of 
symbols per group y = n 2 + g — 1. The weight matrices of 



each Do e 



^nmxnm 



The first block Di 



where 



is 

n X n 



one of the n 
is Hermitian, 
, n are of the 



matrices of the form V ® Ui, where V S 
and the remaining g — 1 matrices Dj, j = 2. 
form ±I n (g) Vj for some set of g unitary m x m matrices 
Ui,...,U s . 

Using an argument similar to the one used with delay- 
optimal codes, it can be shown that the sphere decoding 
complexity of the non delay optimal codes is of the order of 

q V V / / and that the codes are non-singular only 

for M > 2^2- ~ 2i2. For g = 2 and equal values of N and 
T, the non-delay optimal codes of [14| and the codes of |[T2ll 
have the same rate and sphere decoding complexity. 

IV. Discussion 

In this paper we have introduced the notion of singularity 
of STBCs and showed that all known families of high rate 
multigroup ML decodable codes are singular for certain num- 
ber of receive antennas. The following facts which were not 
known before have been shown. 

. Though the N = 4, T = 4 code of [ 10| and the N = 4, 
T = 8 code of [12| have identical rate of iZ cspcu, the 
sphere decoding complexity of the code from ifTUl is less 
than that of the code from lfl2l . 

• For g — 2, and equal values of N, the delay-optimal 
codes of [14 1 and the codes of lfl3l have equal rate and 
the same order of sphere decoding complexity. 

• For equal values of N and T, the codes in Ifl2l and the 
non-delay optimal codes of [14| have identical rate and 
sphere decoding complexities. 

The results and ideas presented in this paper have brought 
to light the following important open problems. 

• Is there an algebraic criterion that ensures that a code 
is non singular? For example, is every code with non 
vanishing determinant also non-singular for all M > RI 

• Do there exist high rate multigroup ML decodable codes 
that are non-singular for arbitrary values of Ml 

• Do there exist singular high rate multigroup ML decod- 
able STBCs with lower sphere decoding complexity than 
that of the known codes? 
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Appendix A 
Proof of TheoremQ] 

From Proposition Q] it is clear that the theorem is true for 
M > N. Thus, we will only consider the case M < N. Before 
giving the proof of Theorem [1] we present two results which 
are used in the proof. Let ei , . . . , ejy be the N columns of the 
matrix Iy, 

Proposition 5: For any i = 1, . . ., N, with probability 1 
(w.p.l), the vector does not belong to the column space 
of the channel matrix H. 

prove the result for i = 1. 
[Hi 
H 2 
Now. 



Proof: 
Let the 

Hi g C N-MXM 

the matrix H 



We first 
channel 



and 



realization H 

H, £ rrMxM 



e c 



NxN 



where 
consider 
We have, 



Ijv-m Hi 
H 2 

det(H) = det(lN_M ) ■ det(H.2) = det(H2) which is non- 
zero w.p.l lfT31 . Thus, w.p.l, the columns of H are linearly 
independent over C, and since ei is the first column of H, 
this means that w.p.l, ei does not belong to the column space 
of the matrix H. 

Now consider any i G {1, . . . , N}. Let P be an N x N 
permutation matrix such that Pe; = ei. Since P is full-ranked, 



Si belongs to the column space of H if and only if ei belongs 
to the column space of PH. Since P is unitary, the distribution 
of H and PH are one and the same, and hence the probability 
that ei belongs to the column space of PH is 0. Thus, with 
probability 1, does not belong to the column space of H. 

■ 

For a given channel realization H, let 
5 = {zeC JV |z H H = 0}. Since H is of rank M w.p.l, 
the dimension of S over R is equal to 2{N — M) w.p.l. Let 
z = [zi Z2 ■ ■ ■ zn] t and for i = 1, . . . , N, let (fi : S — > R 
be the vector space homomorphism that sends the vector z to 
the real number (^)q. We are interested in the dimension of 
the subspace of S which is composed of vectors whose i th 
component is purely real, i.e., in the dimension of ker(ipi). 

Proposition 6: For i = 1, . . . , N, the dimension of image of 
S under the map ipi, dim (ker((fi)) = 2(N — M) — 1 w.p.l. 

Proof: For any given H, fi{S) C R and hence 
dim (<pi(S)) is either or 1. Suppose, dim (tpi(S)) — 0, then, 
there is no vector z in S such that z, is non-zero because 
if such a z exists, the vector iz* ■ z belongs to S and the 
imaginary part of its i th component is \zi\ 2 ^ 0, and thus 
dim((fii(S)) = 1, which is a contradiction. Since, for all the 
vectors in S, the i th component is 0, we have 

S = {z\z H H = 0} = {z\z H [H e,] = 0}. 

Thus, the dimension of the column space of H and [H e^] 
are the same. This means that belongs to the column 
space of H. From Proposition [5] belongs to the column 
space of H w.p.O and hence dim(ipi{S)) = w.p.0. 



Thus, dim((fi(S)) = 1 w.p.l. From rank-nullity theorem, 
dim (ker(ipi)) = dim{S) - dim(cpi(S)) = 2{N - M) - 1 
w.p.l. ■ 

Proof of Theorem Q} Let the weight matrices of the STBC 
QHerm be Ai , . . . , A W 2 and let the space of N x N Hermitian 
matrices over R be given by U = (Ai, . . . , A N 2). For a given 
channel realization H, let p :U — > C JVxM be the R-vector 
space homomorphism that sends the matrix A to AH. Clearly, 
ranfc(G) is equal to the dimension of the subspace p{U) over 
BL Since p(U) is isomorphic to U/ker(p) as vector spaces, 
we have 

ranfc(G) = dim(lA) — dim(ker(p)) = N 2 — dim (ker(p)) . 

Thus, it is enough to show that dim (ker(p)) = (N — M) 2 
w.p.l. Let A £ ker(p) and let a^, . . . , denote the N rows 
of A. Then, ai satisfies af^H = 0, and since the A is Hermi- 
tian, the first component of ai is purely real. From Proposi- 
tion [6] ai £ ker((fi) whose dimension is 2(N —M) — l w.p.l. 
Given a choice of ai, since A is Hermitian, the first component 
of &2 equals the conjugate of the second component of ai, 
and the second component of &2 is purely real. As a result 
of these restrictions and since a^H = 0, i.e., &2 € ker((fi2), 
&2 belongs to a coset of the subspace of S whose dimension 
is (2(N - M) - 1) - 2 = 2{N - M) - 3. Similarly, given a 
choice for ai, . . . , &k-i, where k = 1, . . . , N, the first k — 1 
components of a& are fixed and the imaginary part of the 
k th component is zero. Hence, a^ belongs to a coset of a 
subspace of S with dimension (2(N — M) — 1) — 2k. Thus, 
the dimension of ker (p) equals 

2(N-M)-1 + 2{N-M)-3 + ••• + 1 = (N-M) 2 , 
with probability 1. This completes the proof. 
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